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Abstract: The origin of -R-parity in supersymmetric models can be explained if 
B — L is part of the gauge group. We discuss the mass spectrum of the minimal 
U(l) Y x U(1) B -l model based on a GUT implementation using CMSSM-like bound- 
ary conditions. Here we focus in particular on the Higgs and neutralino sectors in 
this class of models. While the neutralinos can have masses as low as 100 GeV, we 
show that the requirement of being consistent with existing bounds on the Z' implies 
that in general the sfermions have masses in the multi-TeV range. In the extended 
Higgs sector we show the existence of a second light state which, however, will be 
difficult to observe, while having at the same time a SM-like Higgs in a mass range of 
123-126 GeV. Moreover, we propose a set of benchmark scenarios for phenomenolog- 
ical studies. On the technical side we demonstrate that gauge kinetic mixing effects 
can be quite important, affecting in particular the Higgs and the neutralino sectors. 
Not only can they shift the mass of the lightest neutralino by about 10 per-cent but 
also they can change the nature of neutralinos and Higgs bosons in a significant way. 
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1. Introduction 

Models with an additional U{1)b~l gauge symmetry at the TeV scale have recently 
received considerable attention. On one hand, they are among the simplest extensions 
of Standard Model (SM) gauge group with observable consequences at the LHC 
[1, 2, 3, 4]. On the other, this class of models can help to understand the origin of 
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.R-parity and its possible spontaneous violation in supersymmetric models [5, 6, 7], 
as well as the mechanism of leptogenesis [8, 9]. It has been shown that a gauge sector 
containing U(1)y x U(1)b-l can be a result of an E$ x E 8 heterotic string theory 
(and hence M-theory) [10]. While most studies of supersymmetric variants have so 
far focused on the effects of the additional gauge group far below the GUT scale, the 
questions arise of whether this group can be unified at the high scale with the SM 
gauge group and what the phenomenological consequences are. A renormalization 
group equation (RGE) analysis of such a model, assuming the unification of the 
gauge groups, has been performed in [7, 11]. However, the effects of possible mixing 
between the two Abelian groups have been neglected so far: it is well known that in 
models with several U(l) gauge groups, kinetic mixing terms 

-XabF^fy, a^b (1.1) 

between the field strength tensors are allowed by gauge and Lorentz invariance [12], 
as F a ^ v and F h ^ u are gauge invariant quantities by themselves, see e.g. [13]. Even 
if these terms are absent at tree level at a particular scale, they might be generated 
by RGE effects [14, 15]. 

The impact of gauge kinetic mixing in generic extensions of the standard model 
(SM) and the MSSM has been studied so far in several aspects. For instance, one 
can show that the dark matter of the universe can be charged with respect to an 
additional U(l) but neutral with respect to the SM gauge group. However, here one 
can show that there is a residual SM gauge interaction of the dark matter particles 
due to the gauge kinetic mixing. The consequences for the relic density and the cross 
sections concerning direct as well indirect detection of dark matter have been ana- 
lyzed [16, 17, 18]. It has been shown that these cosmological bounds are sometimes 
more severe than the bounds from electroweak precision data if the dark matter can- 
didate interacts dominantly due to kinetic mixing. Moreover, the kinetic mixing in 
the context of supersymmetric hidden sector dark matter has been considered in [19] 
and the LHC phenomenology of a nearly decoupled sector only interacting with the 
visible sector due to kinetic mixing has been elaborated in [20] . 

In this work, we discuss the mass spectrum of the model presented in [5, 7]. This 
minimal B — L extension of the Minimal Supersymmetric Standard Model (MSSM) 
has a U{\)b-l gauge group tensored to the SM gauge groups and two bileptonic 
chiral superfields which are gauge singlets under SM gauge groups. In addition, 
three right-handed neutrinos are needed to ensure that U(1)b-l is anomaly-free, 
which provide the necessary ingredients to explain neutrino data. We refer to this 
model as the BLSSM. 

The focus of this paper is on the mass spectrum of this model and resulting 
phenomenological aspects assuming mSUGRA-like boundary conditions at the GUT 
scale and unification of the B — L coupling with the SM couplings. In particular we 
will demonstrate that gauge kinetic mixing effects are particularly important in the 



2 



Higgs and neutralino sectors. These effects do not only change the masses of these 
particles but have quite some impact on their nature, e.g. they induce tree-level 
mixing which would be absent if these effects were to be neglected. We will show 
that new light Higgs states are possible without being in conflict with current data 
while having at the same time a SM-like Higgs in the range close to 120 GeV. We 
will focus here on the case of R-parity conservation and discuss the case of broken 
R-parity violation in a subsequent paper [21]. 

In the usual CMSSM with the MSSM particle content, the lightest neutralino 
is mainly bino-like. We show that in our model the nature of this particle can 
be quite different and identify regions where it is either mainly a S , f/(2)x,-doublet 
Higgsino, a £/(l)£_£-gaugino which we dub the BLino, or a fermionic partner of the 
^(l)B-L-breaking scalar which we dub the bileptino, since we call the scalar the 
bilepton for reasons given below. In the next section we introduce the model and 
focus in particular on aspects related to the spectrum. In section 3 we present our 
numerical results and provide benchmark points with distinct features and in section 
4 we draw our conclusions. In the appendices we collect supplementary formulas 
for mass matrices, anomalous dimensions and /3-functions at lowest order needed for 
the discussion of the main features in section. The corresponding formulas including 
higher order effects can be easily computed using the input files for SARAH given in 
appendix C. 

2. The Model 

In this section we present the particle content of the model considered. An important 
aspect is the U(l) gauge kinetic mixing which is discussed in some detail as it leads 
to significant changes in the spectrum. Although we include loop corrections for the 
numerical analysis when calculating the masses, we restrict ourselves in this section 
to tree-level expressions, as this is sufficient for discussing the main differences with 
respect to the MSSM. 

2.1 Particle content and superpotential 

The model consists of three generations of matter particles including right-handed 
neutrinos which can, for example, be embedded in 5*0(10) 16-plets. Moreover, below 
the GUT scale the usual MSSM Higgs doublets are present as well as two fields i] 
and fj responsible for the breaking of the U(1)b-l- Furthermore, rj is responsible 
for generating a Majorana mass term for the right-handed neutrinos and thus we 
interpret the B — L charge of this field as its lepton number, and likewise for f/, and 
call these fields bileptons since they carry twice the lepton number of (anti-) neutrinos. 
We summarize the quantum numbers of the chiral superfields with respect to U (1) y x 
SU(2) L x SU(3) C x U(1) B -l in Table 1. 
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L 


L 


3 


(-- 2 1 -i) 

v 2 7 7 7 2 7 
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e c 


3 


(1,1,1,1) 


V 


u c 


u c 


3 


(0,1,1,|) 


H d 


H d 


H d 


1 


(-1,2,1,0) 


H u 


H u 


H u 


1 


(1,2,1,0) 


fj 


V 


fj 


1 


(0,1,1,-1) 


fj 


fj 


fj 


1 


(0,1,1,1) 



Table 1: Chiral superfields and their quantum numbers. 



The superpotential is given by 

W=Yi*UiQjH u - Y^ Di Qj H d -Y^E t L 3 H d + fiH u H d 

+ Y l J U H u uj -fx'fjfj + Y$ h f] u 3 (2.1) 

and we have the additional soft SUSY-breaking terms: 

£sb =£mssm - \b^b> m bb> - ^\ b ,\b,M b > - m 2 v \r]\ 2 - m\\f)^ - m 2 vi ^{vffv] 

- rjfjB^ + TyHJZLj + y;V' "- (2.2) 

i,j are generation indices. Without loss of generality one can take and to be 
real. The extended gauge group breaks to SU(3) C <8> U(l) em as the Higgs fields and 
bileptons receive vacuum expectation values (vevs): 

H ° d= 7^> ^ + Vd + l<t)d) ' H u= ^^u + v u + i<p u ) (2.3) 
r) = -^=(a v + v ri + i(f) v ) , fj = fa + Vfj + ifc) (2.4) 

We define tan/3' = jjjj in analogy to the ratio of the MSSM vevs (tan/3 = ^). 
2.2 Gauge kinetic mixing 

As already mentioned in the introduction, the presence of two Abelian gauge groups 
in combination with the given particle content gives rise to a new effect absent in the 
MSSM or other SUSY models with just one Abelian gauge group: the gauge kinetic 
mixing. This can be seen most easily by inspecting the matrix of the anomalous 
dimension, which at one loop is given by 
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where the indices a and b run over all U(l) groups and the trace runs over all fields 
charged under the corresponding U(l) group. 
For our model we obtain 

7=^("fU- ( 2 - 6 ) 



16tt 2 V 4 6 



and we see that there are sizable off-diagonal elements. N contains the GUT nor- 
malization of the two Abelian gauge groups. We will take as in ref. [7] for U(l)y 

and J\ for U(1)b-l, i-e. N = diag(w|, y §)■ Hence, we obtain finally 




7 = is? ' - ^ ' ' < 27) 



Therefore, even if at the GUT scale the U(l) kinetic mixing terms are zero, they are 
induced via RGE evaluation at lower scales. In practice it turns out that it is easier 
to work with non-canonical covariant derivatives instead of off-diagonal field-strength 
tensors such as in eq. (1.1). However, both approaches are equivalent [22]. Hence in 
the following, we consider covariant derivatives of the form 

D li = d ll -iQlGA (2.8) 

where is a vector containing the charges of the field <fi with respect to the two 
Abelian gauge groups, G is the gauge coupling matrix 

9yy 9yb \ 

9by 9bb . 



G= [ » rr (2.9) 
) 



and A contains the gauge bosons A = (A?, A^) T . 

As long as the two Abelian gauge groups are unbroken, we have still the freedom 
to perform a change of basis: A = (A^,A^) ->• A = ((A^)', (A^)') = RA where 
R is an orthogonal matrix. It is possible to absorb this rotation of the gauge fields 
completely in the definition of the gauge couplings without the necessity of changing 
the charges, which can easily be seen using eq. (2.8) 

QlGA = QlG(R T R)A = Q^(GR T )A' = Q%GA' (2.10) 

This freedom can be used to choose a basis such that electroweak precision data can 
be accommodated in an easy way. A convenient choice is the basis where qby = 
as in this basis only the Higgs doublets contribute to the entries in the gauge boson 
mass matrix of the U(l)y <8> SU(2)l sector and the impact of r\ and fj is only in the 
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off-diagonal elements as discussed in section 2.4. Therefore we choose the following 
basis at the electroweak scale [23]: 

/ 9yy9bb — 9yb9by /<-, i i \ 

9yy= / 2 = 2 ~ = 9i (2.11) 

V 9bb + 9by 



9'bb=\J9bb + 9 by = 9bl (2.12) 
/ 9yb9bb + 9by9yy ~ /oiq\ 

9YB = / 2 , 2 ' = 9 ^- 13 ) 

V 9bb + 9 BY 

9by =0 (2.14) 
This also leads to our condition for finding the GUT scale in the numerical analysis: 

9yy9bb — 9yb9by /01 ^ 

92 = /2 ^ 2 ( 2 - 15 ) 

V 9bb + 9 BY 

This is equivalent to a rotation of the general 2x2 gauge coupling matrix at each 
energy scale to the triangle form and using g\ = g 2 as the GUT condition. Neglecting 
threshold corrections, this leads in the case of kinetic mixing to exactly the same GUT 
scale as in the MSSM [24]. 

Immediate interesting consequences of the gauge kinetic mixing arise in various 
sectors of the model as discussed in the subsequent sections: (i) it induces mixing at 
tree level between the H u , H d and rj, fj; (ii) additional D-terms contribute to the mass 
matrices of the squarks and sleptons; (iii) off-diagonal soft-SUSY breaking terms for 
the gauginos are induced via RGE evolution [22, 25] with important consequences for 
the neutralino sector as discussed in section 3.4, even if at some fixed scale M ab = 
for a b. 

2.3 Tadpole equations 

We find for the four minimization conditions at tree level 

td = v d (m 2 Hd + l/f + I (g\ + gl + g 2 ) (v\ - v 2 u ) + ±gg BL (v 2 v -v 2 )^- v u B^ = 

(2.16) 

tu = v u [m 2 Hu + |/f + ^ {g\ + g\ + g 2 ) (v 2 u - v 2 d ) + ^gg B L (uj - v 2 )^j - v d B^ = 

(2.17) 

t v = v v (m 2 + |/if + ±gg BL (v 2 d - v 2 u ) + ^g 2 BL (v 2 - v 2 ^ - v n B^ = (2.18) 
h = v n (m 2 + |/if + ±gg BL (v 2 u - v 2 d ) + X -g 2 BL (v 2 - v 2 )^j - v v B„, = (2.19) 

We solve them with respect to /i, -B M ,// and B^ as these parameters do not enter 
any of the RGEs of the other parameters. Using x 2 = v 2 + t>| and v 2 = v 2 d + v 2 u we 
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obtain 

|/x| 2 =g [{2~gg BL x 2 cos(2/3') - 4m^ + 4m|J sec(2/3) - ±(m 2 Hi + to 2 ^) - (<? 2 + ~ g 2 + 

(2.20) 

S f = - ^ ( - ^99bl^ cos(2/3') + 4m^ d - 4m^ + ^ + f + cos ( 2/3 ) j t an(2/3) 

(2.21) 

|/if =1 ( - 2 ^| L x 2 + m 2 + m 2 ) + (2m 2 - 2m 2 + qqbl^ cos(2/3)) 860(2/3')) 

(2.22) 

=- ( - 2g 2 BL x 2 008(2/3') + 1m\ - 2m 2 + qqbl^ cos(2/3)) tan(2/3') (2.23) 

~ Qblx as we will show in section 2.4 and, thus, we find an approximate relation 
between M' z and // 

4(m 2 -m 2 tan 2 /3') - t; 2 ^ BL cos/3(l + tan/3') 
" - 2|/i 1 + 2(tan 2 /3'-l) (2 ' 24) 

A closer inspection of the system shows that either m? or m 2 has to become negative 
to break U (1)b-l- For both parameters, gauge couplings enter the RGEs, increasing 
their values when evolving from the GUT scale to the electroweak scale. The Yukawa 
couplings Y v and Y x as well as the trilinear couplings T u and T x lead to a decrease, 
but at the one-loop level they only affect the RGE for m 2 . However, neutrino data 
require \Y u>i j\ to be very small in this model and thus they can be neglected for 
these considerations. Therefore, will always be positive whereas m 2 can become 
negative for sufficient large Y x and T x . 

We can roughly estimate the contribution of these couplings to the running 
value of m 2 by a one-step integration assuming mSUGRA-like GUT conditions (see 
sect. 2.8) to 

^ = ~-a(Y^)(3ml + 2© log (2.25) 

with T x ~ A Y X . Therefore, we expect that large values of m and A will be 
preferred, implying heavy sfermions. Moreover, tan/3' has to be small and of 0(1) in 
order to get a small denominator in the second term of eq. 2.24. One last comment 
concerning the effect of gauge kinetic mixing: qyb is always negative below the GUT 
scale if it is zero at the GUT scale as can be seen by the following: for vanishing 
off-diagonal gauge couplings, the /3-functions eqs. (B.17) and (B.16) will always be 
positive, i.e. qby and qyb are driven negative. Using eq. (2.13), one can see that this 
also drives g negative. Therefore, the second term will give a positive contribution. 
From this point of view, one might expect that small m for given tan/3' would be 
sufficient to get the same size of However, as can been seen in Fig. 1, where we 
plot the tree-level value of // in the (m ,tan/3')-plane for the cases with and without 
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m [GeV] m [GeV] 



Figure 1: Contour plots of [/,' at tree-level in the (mo, tan/3')-plane for Mz> = 2000 GeV 
(left) and M z > = 4000 GeV (right). The other parameters are M 1/2 = 0.5 TeV, tan(/3) = 10, 
Aq = 1.5 TeV, Y Xj u = 0.42. The full lines correspond to the case including gauge kinetic 
mixing, the dashed lines are without kinetic mixing. 

kinetic mixing, the opposite effect takes place. The reason is the contribution of the 
kinetic mixing to the evaluation of and m?. One can also see in this figure that 
the upper limit of tan/3' for a given value of mo decreases with increasing Mz> as 
expected. Even if one might get the impression from this figure that the effects of 
kinetic mixing are in general small as they slightly shift the region where breaking 
of U(1)b-l can occur, it will be shown later that it can have a significant impact on 
the masses. 

For the numerical results we include one-loop corrections to eqs. (2.16)-(2.23) 
as well as for all masses. This is done by using the DR scheme and extending the 
MSSM results given in ref. [26] in a similar manner to the NMSSM case discussed in 
ref. [27]. We denote the one-loop contributions to the tadpole equations (2.16)-(2.19) 
by St!- . The requirement of keeping the values of tan/3 and tan/3' after including 
the loop corrections as well as the conditions 

ti + 5t[ 1] = for % = d, u, T], fj (2.26) 

lead to shifts of /i, B^, // and B^i compared to the values obtained by eqs. (2.20)- 
(2.23). 

2.4 Gauge boson mixing 

Due to the presence of the kinetic mixing terms, the B' boson mixes at tree level 
with the B and W 3 bosons. Requiring the conditions of eqs. (2.11)-(2.14) means 
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that the corresponding mass matrix reads, in the basis (B, W 3 , _£>'), 









l-( 


\b>) 





\gW -\g1g2v 2 \gi~gv 2 
-\g1g2v 2 \g\v 2 -\~gg 2 v 2 | (2.27) 
\gigv 2 -\~gg2v 2 (g 2 BL x 2 + \g 2 v 2 ) 

In the limit g — > both sectors decouple and the upper 2x2 block is just the 
standard mass matrix of the neutral gauge bosons in EWSB. This mass matrix can 
be diagonalized by a unitary mixing matrix to get the physical mass eigenstates 7, 
Z and Z' . The rotation matrix can be expressed by two mixing angles ®w and Q' w 

as 

cos ®w cos Q'w sin 6w — sin Qw sin Q'w \ / 7 \ 
sin Q w — cos Qw cos Q'w cos0jysin0V ^ (2.28) 

sinGV cosGV / \Z' J 

The third angle is zero due to the special form of this matrix. Q'w can be approxi- 
mated by [28] 

tan 20V * 2 _tM±M (2.29) 

+ 16(f) \g 2 BL -gl-g\ 

The exact eigenvalues of eq. (2.27) are given by 

M 7 = (2.30) 
v = \ ((<7i 2 + ^2 + fV + 4g 2 BL x 2 T 

+ si + ^ 2 ) 2 ^ 4 - %? + g 2 2 )g 2 BL v 2 x* + 16^) (2.31) 

Expanding these formulas in powers of v 2 /x 2 , we find up to first order: 

Mz = \ {gl + gl) v 2 , M z , = g 2 BL x 2 + (2.32) 

All parameters in eqs. (2. 16)- (2. 19) as well as in the following mass matrices are 
understood as running parameters at a given renormalization scale Q. Note that the 
vevs Vd and v u are obtained from the running mass Mz(Q) of the Z boson, which is 
related to the pole mass M z through 

M 2 Z (Q) = 9 -^(v 2 u + v 2 d ) = M 2 z + Re{U T zz (M 2 z )}. (2.33) 

Here, Il^ z is the transverse self-energy of the Z. See for more details also ref. [26]. 

The mass of additional vector bosons as well as their mixing with the SM Z 
boson, which imply for example a deviation of the fermion couplings to the Z boson 
compared to SM expectations, is severely constrainted by precision measurements 



M 



z,z> 
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from the LEP experiments [29, 30, 31]. The bounds are on both the mass of the 
Z 1 and the mixing with the standard Z boson, where the latter is constrained by 
| sm(Qw') < 0.0002 1 . Using eq. (2.29) together with eq. (2.32) as well as the values 
of the running gauge couplings, a limit on the Z' mass of about 1.2 TeV is obtained. 
Taking in addition the bounds obtained from U, T and S parameters into account 
[32] one gets — > 7.1 TeV which for g Bh ~ 0.52 implies M z > > 1.8 TeV. 

Ve QB—L 

Therefore we have taken always M z > > 2 TeV. In this way we have also satisfied the 
most recent bounds obtained by the ATLAS and CMS [33]. 



2.5 The Higgs sector 

In this section we present the tree-level formulas for the Higgs sector and we briefly 
discuss the main steps to include the one-loop corrections. The one-loop formulas 
and further details will be presented elsewhere [34]. 



2.5.1 Pseudoscalar Higgs bosons 

It turns out that in this sector there is no mixing between the 577(2) doublets and 
the bileptons at tree level and we obtain in the basis (4>d, <f> u , <f>r), <f>fj)' 



I B u tan 8 B 



m 



A,T 



\ 










B^ cot 8 
Bp. tan 8' B„ 



\ 







B, 



Bp, cot 8' J 



(2.34) 



Obviously, both sectors decouple at tree level. This is a consequence of the fact that 
we assume that there is no CP violation in the Higgs sector, so the different D-term 
contributions cancel exactly. One obtains two physical states A and A® with masses 



m 2 A o 



2£ M 
sin 28 ' 



2 2iV 



(2.35) 



A more detailed study of the pseudoscalar sector at one-loop, including the question 
if the block-diagonal form of the mass matrix in eq. (2.34) can be maintained at 
higher order, goes beyond the scope of this work and will be presented elsewhere 
[34]. 



2.5.2 Scalar Higgs bosons 

In the scalar sector the gauge kinetic terms do induce a mixing between the SU (2) 
doublet Higgs fields and the bileptons. The mass matrix reads at tree level in the 
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basis (cr d ,(7 u ,(7^,(T^): 



rn 



h,T — 

-m 2 A0 cpSp - g 2 v d v u 

_ 99BL r. 



-m^cpsp 
m 2 A0 c 2 + 

99 BL 



- g 2 VdVv 
g 2v d 



■-V d Vfj 



2 

99BL 
2 



V u Vri 



V u Vfj 



^VdV, 

99 BL „, „, 
2 u u u r) 

7TI 'aoc 2 / + g 2 BL v % 



99BL 



VdVfj 



'-11 11- 



-mjocpsp 



gBL v n v v 



2 

99BL , 
2 



where we have defined g 2 = \ (g 2 + g\ + g 



7,2\ 



cosfa;) and s T = sinfx) 



V 



(2.36) 



The one-loop corrections are included by calculating the real part of the poles of the 
corresponding propagator matrices [26, 34] 



Bet [p 2 l - ml lL (p 2 )} = 0, 



where 



m 



h,lL 



(P 2 



2,h 



n hh (p 2 



Equation (2.37) has to be solved for each eigenvalue p 2 
in an iterative procedure. 



(2.37) 
(2.38) 

m 2 which can be achieved 



2.5.3 The charged Higgs boson 

At the tree level one finds that the charged Higgs boson mass has exactly the same 
form as in the MSSM: 

m 2 H+ = 171% + mly (2.39) 

However, for the one-loop corrections one obtains additional contributions due to the 
kinetic gauge mixing [34]. 

2.6 Neutralinos 

In the neutralino sector we find that the gauge kinetic effects lead to a mixing between 
the usual MSSM neutralinos with the additional states, similar to the mixing in the 
CP-even Higgs sector. In other words, were these to be neglected, both sectors would 
decouple. The mass matrix reads in the basis ^A^, W°, H®, H®, \ B ,, fj, ijj 



/ 


Mi 





-\g\Vd 


\g\v u 


\M BB , 





\ 







M 2 


\g2Vd 


-\g2Vu 












-\g\Vd 


\g2Vd 







-\gv d 









\g\v u 


-\g2Vu 


-n 





\gv u 









\M BB , 





-\gvd 


\gv u 


M B 


-gBLV v gBLVfj 


















-gBLV v 


-fi' 




\ 














gBLVfj 


-// 


/ 



(2.40) 
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It is well known that for real parameters such a matrix can be diagonalized by an or- 
thogonal mixing matrix N such that N*M£ W is diagonal. For complex parameters 
one has to diagonalize M£ (M£ y . We obtain, in a straightforward generalization 
of the formulas given in [26], at the one-loop level 



ES(P?) + E°/(P?) + (e°' T (p?) + ^(P?)) M|° 
+M|° (ES T (p?) + E£(p?)) 



(2.41) 



where we have denoted the wave-function corrections by Y,° R , TP L and the direct one- 
loop contribution to the mass by TP S . 

In this model, for the chosen boundary conditions, the lightest supersymmetric 
particle (LSP), and therefore the dark matter candidate, is always either the lightest 
neutralino or the lightest sneutrino. The reason is that m must be very heavy in 
order to solve the tadpole equations, and therefore all sfermions are heavier than the 
lightest neutralino, with the possible exception of the sneutrinos. As described in 
sec. 2.7 below, the splitting of the CP-even and CP-odd components of the sneutrinos 
can be very large, pushing the mass of the lighter eigenstate down even to the point 
of being lighter than the lightest neutralino. However, we have chosen to leave 
the investigation of such a scenario to future work, and we take only points with 
neutralino LSPs as our benchmark scenarios. Before leaving this topic, we note 
that BLV has a lightest sneutrino that is almost as light as the lightest neutralino. 
A neutralino LSP is in general a mixture of all seven gauge eigenstates. However, 
normally the character is dominated by only one or two constituents. In that context, 
we can distinguish the following extreme cases: 

1. Mi < M 2 ,/jl,M b ,ijl': Bino-like LSP 

2. M 2 < M x ,n,M B ,n': Wino-like LSP 

3. [i < M 1 ,M 2 ,M B ,/ji': Higgsino-like LSP 

4. M B < Mi,M 2 ,/x,//: BLino-like LSP 

5. jj! <C Mi, M 2 , /i, M B : Bileptino-like LSP 

In addition, we will summarize the Bino- and Wino-like states, i.e. the states built 
by the gauginos of the MSSM, in the following often by 'gaugino-like'. Note that 
this doesn't include the BLino, the gaugino of the B — L sector. 

Although the gauge kinetic effects do lead to sizable effects in the spectrum, they 
are not large enough to lead to a large mixing between the usual MSSM-like states 
and the new ones. Therefore, we find that the LSP is either mainly a MSSM-like 
state or mainly an admixture between the BLino and the bileptinos. A discussion of 
the parameter space where the different characters appear is given in sec. 3.4. 



12 



2.7 Charginos and sfermions 

For completeness we also give a short summary of the other sectors of the model. 
The chargino mass matrix at tree level is exactly the same as for the MSSM: 

Mf=(\ M2 * 92 ""V (2.42) 

This mass matrix is diagonalized by a biunitary transformation such that U*M^, 
is diagonal. The matrices U and V are obtained by diagonalizing M% (M£ )t and 
(Mj, )*{Mj, ) T , respectively. At the one-loop level, one has to add the self-energies 
[26] 

m£( P >) = Mf - E+(rf) - E+(p?)JWf - M^tipl). (2.43) 

The mass matrices for the squarks and the charged sleptons are given in ap- 
pendix A. At tree level, the differences compared to the MSSM are the additional 
D-terms in the diagonal entries. All complex scalar mass matrices are diagonalized 
by an unitary matrix Z 

Z+m\z\ = ml diag . (2.44) 

The corresponding mass matrices at the one-loop level are again obtained by taking 
into account the self-energy according to 

mlt(p*)=m 2 /-Tl H (pt), (2.45) 

and the one-loop masses are obtained by calculating the poles of the real part of the 
propagator matrix. 

We focus here on the sneutrino sector as it shows two distinct features compared 
to the MSSM. Firstly, it gets enlarged by the additional superpartners of the right- 
handed neutrinos. Secondly, even more drastically, a splitting between the real and 
imaginary parts of the sneutrino occurs resulting in twelve states: six scalar sneutri- 
nos and six pseudoscalar ones [35, 36]. The origin of this splitting is the Y£* Vif]Vj 
in the superpotential, eq. (2.1), which is a AL = 2 operator. Therefore, we define 

i>l = ^= K + Wl) ?r = ^= K + (2-46) 

The 6x6 mass matrices of the CP-even and CP-odd sneutrinos can be written in 
the basis (o\l,or) respectively (4>l,<Pr) as 

ml^=x( mLLm * L ) (2.47) 
V m RL m RR I 
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with 

mLL = \( 1 ( (al + gl + g 2 )(-vl + vfj + ~gg B (- H, + H - v 2 u + vfj 



+ 2g 2 B (-v 2 + v 2 ) ) + 8m? + toftrtf) (2-48) 
m*L = l - ( - 2v^u d ^ + v u (2V2TJ ± ^^y;) ) (2.49) 
m RR = l -(l (2g 2 B (v 2 n - v 2 ) -gg B (-v 2 u + vfj ) + 8ml + 2 Vfj ( =F ±y/2Y x n'*) 

+ Av 2 u Yj y; + 2^ ( ± 4^ + 8^7 x j) (2.50) 



where we have assumed CP conservation. In the case of complex trilinear couplings, 
a mixing between the scalar and pseudoscalar particles occurs, resulting in 12 mixed 
states and consequently in a 12 x 12 mass matrix. It particular the term ~ VfjY x ji'* 
is potentially large and induces a large mass splitting between the scalar and pseu- 
doscalar states. Also the corresponding soft SUSY term ~ v v T x can lead to a sizable 
mass splitting in the case of large \Aq\. As a side-remark we note that in such a case 
also the determinant of this matrix could become negative indicating the breaking 
of R-parity in a somewhat different way compared to the discussion in ref . [7] . How- 
ever, here we will concentrate on the R-parity conserving case and discuss R-parity 
violation elsewhere [21]. 

2.8 Boundary conditions at the GUT scale 

We will consider in the following a scenario motivated by minimal supergravity. This 
means that we assume a GUT unification of all soft-breaking scalar masses as well 
as a unification of all gaugino mass parameters 

22222 /n ri\ 

m o = m H d = m H u =m v = m n (2.51) 
mlSij =m 2 D 5ij = mljSij = m^y = m 2 E 8ij = m 2 L 5ij = m 2 v 8ij (2.52) 
Mi/a =Mi = M 2 = M 3 = M B , (2.53) 

Also, for the trilinear soft-breaking coupling, the ordinary mSUGRA conditions are 
assumed 

Ti = A Yi, i = e,d,u,x,v . ( 2 -54) 

We do not fix the parameters /x, // and B^i at the GUT scale but determine 
them from the tadpole equations. The reason is that they do not enter the RGEs 
of the other parameters and thus can be treated independently. The corresponding 
formulas are given in section 2.3. 

In addition, we consider the mass of the Z' and tan/3' as inputs and use the 
following set of free parameters 

mo, M1/2, Aq, tan/3, tan/3', sign(/i), sign(//), M z >, Y x and Y v . (2.55) 
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Figure 2: One- loop evaluation of the gauge couplings (left) and gaugino mass parameters 
(right) associated with the Abelian gauge groups. Blue color is used for gyy & Mbb, black 
for gsB & Mb'b'i dashed black for gys & Mbb 1 and dashed blue for gsY- The evolution 
of the two off-diagonal couplings hardly differs and thus the two lines nearly match. 



Y v is constrained by neutrino data and must therefore be very small in comparison 
to the other couplings; thus it can be neglected in the following. Y x can always be 
taken diagonal and thus effectively we have 9 free parameters and two signs. 

Furthermore, we assume that there are no off-diagonal gauge couplings or gaug- 
ino mass parameters present at the GUT scale 



9by =9yb = 
M BB > =0 



(2.56) 
(2.57) 



This choice is motivated by the possibility that the two Abelian groups are a remnant 
of a larger product group which gets broken at the GUT scale as stated in the 
introduction. In that case gyy and qbb correspond to the physical couplings g\ and 
Qbli for which we assume a unification with g 2 : 



g GUT 



g^ T 



9bl 



(2.58) 



where we have already taken into account the correct GUT normalization as dis- 
cussed in section 2.2. In Fig. 2 we display the running of the gauge couplings and 
gaugino parameters in the Abelian sector to demonstrate the effect of the gauge 
kinetic mixing. The GUT scale has been set to 2 • 10 16 GeV where we have fixed 
9yy = 9bb = 0.72 and Mbb = Mb>b' = 200 GeV. All off-diagonal parameters have 
been set to zero. Note that in particular Mbb 1 becomes sizable at the electroweak 
scale. 



3. Numerical results 

All analytic expressions for masses, vertices, RGEs as well as one-loop corrections to 
the masses and tadpoles were calculated using the SARAH package [37, 38, 39]. For 
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the generic expressions, those of ref. [40] are used in the most general form respecting 
the complete flavour structure. In addition, gauge kinetic mixing effects in the RGEs 
are included using the extensions of ref. [22]. The loop corrections to all masses as 
well as to the tadpoles are derived in DR scheme and the 't Hooft gauge. 

The numerical evaluation of the model is very similar to that of the default 
implementation of the MSSM in SPheno [41, 42]: as the starting point, the SM 
gauge and Yukawa couplings are determined using one-loop relations as given in 
ref. [26] which are extended to our model. The vacuum expectation values Vd and 
v u are calculated with respect to the given value of tan/3 at Mz, while v v and are 
derived from the input values of Mz> and tan/3' at the SUSY scale. 

The RGEs for the gauge and Yukawa couplings are evaluated up to the SUSY 
scale, where the input values of Y v and Y x are set. Afterwards, a further evaluation 
of the RGEs up to the GUT scale takes place. After setting the boundary conditions 
all parameters are evaluated back to the SUSY scale. There, the one-loop-corrected 
SUSY masses are calculated using on-shell external momenta. These steps are iter- 
ated until the relative change of all masses between two iterations is below 1CT 4 . 

3.1 Benchmark Points 

For the numerical analysis in the following we have chosen five points in the con- 
strained parameter space which have distinct features. An overview of these points 
is given in Tab. 2. A typical feature is that m has to be in the TeV range to be 
consistent with the existing bounds on the Z'-mass. 

BLI provides a neutralino LSP similar to the MSSM and also the lightest scalar 
Higgs is very similar to that of the MSSM and serves mainly to exemplify technical 
details. BLII demonstrates that it is possible to increase the Higgs mass of the 
MSSM-like light Higgs boson through a mixing with the new B — L fields. BLIII has 
an LSP with a large Higgsino fraction that is difficult to reach within the MSSM. 
BLIV and BLV show two new dark matter candidates: a BLino and a bileptino LSP. 

We have also chosen points which lead to SM-like Higgs masses in the preferred 
range of 123-126 GeV [43, 44]. As we will show in the following, the extended 
Higgs sector has an impact also the the MSSM Higgs masses. For instance, using 
the CMSSM parameters of BLI, the light Higgs would have a mass of 122.2 GeV, 
i.e. 1.5 GeV lighter. However, this effect will be smaller for points like BLIII or BLIV 
for which a larger mass splitting between the Higgs fields of both sectors is present. 
In this case, the light Higgs mass in the BLSSM agrees with that of the MSSM for 
the same parameters. One can also see that the branching ratios to two photons are 
always larger than in the SM except for BLIV, where the additional neutralinos are 
so light that Br(/i — > XiXi) = 43.6%. This feature significantly softens, of course, 
the bounds on the Higgs mass. 
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Figure 3: Spectrum of BLII (upper left figure), BLIII (upper right figure), BLIV (lower 
left figure) and BLV (lower right figure). All masses are given in GeV. The masses of the 
neutralinos and scalar Higgs fields are also given in table 2. 



3.2 Precision of the mass calculation: the impact of kinetic mixing 

Before we take a closer look on the Higgs and neutralino sector, we want to comment 
on the precision of the mass calculation. As already mentioned, we use two-loop 
RGEs and one-loop corrections to the masses. In addition, we take the feature 
of kinetic mixing into account which is often neglected in literature. To show the 
importance of the kinetic mixing and to compare the resulting mass spectrum of 
the B — L model with the MSSM, we show in Tab. 3 the masses calculated at the 
one-loop level using one- or two-loop RGEs. For an easier comparison of the various 
effects, we fix the gauge the GUT scale to 2 • 10 16 GeV and the gauge couplings at 
the GUT scale to g\ = g 2 = g% = Qbl = 0.72. 

It can be seen that the most pronounced differences are in the neutralino and 
Higgs sectors. In case of sfermions and the heavier Higgs bosons the differences 
between the MSSM and BLSSM particle spectrum is about 1-2 percent and thus 
relatively small. They are larger in the charged slepton sector since here the ad- 
ditional U(1)b-l leads to larger effects. The main reason for the smallness of the 
differences is the required largeness of m . Moreover, the D-term effects due to the 
extra U(1)b-l are also small as tan /3' is close to 1. In the sneutrino sector the scalar 



17 




700 800 900 1000 1100 700 800 900 1000 1100 1200 



m [GeV] m [GeV] 

Figure 4: Mass of the lightest Higgs (first row) and the logarithm of the BLino fraction 
of the lightest Higgs (second row). The other parameters are those of BLI. First left: with 
kinetic mixing, right: without kinetic mixing. 

and pseudoscalar particles behave quite differently: while the mass shifts in the pseu- 
doscalar sector are rather moderate, the masses of the scalars change by more than 
100 GeV. The cause for this is a partial cancellation for the given parameter point 
between the large positive terms in eq. (2.50), ml + v v Y x Y*, and the large negative 
term — 4"*/2Y x /j!* which is very sensitive to the exact values of all parameters at the 
SUSY scale. This also implies that kinetic mixing is particularly important for this 
sector and could even trigger i?-parity violation [21]. 

In this table we also demonstrate that neglecting the gauge kinetic effects can 
have dramatic effects in the extended Higgs sector: the mass of hi state would be 
predicted to be about 20% larger. 

Also in the neutralino sector there is quite some impact for the LSP masses 
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which would be wrong by about 10% if gauge kinetic terms were to be neglected. 
But also the properties of the lightest Higgs particle can change: we show in Fig. 4 a 
comparison between the mass and bilepton fraction of the lightest with and without 
kinetic mixing. It can be seen that the masses are clearly shifted while, of course, 
there is also huge difference of several orders in the bilepton fraction between both 
cases. While the bilepton contribution for MSSM-like scalars in the case without 
kinetic mixing is solely based on the mixing at one-loop level, the off-diagonal gauge 
couplings introduce already a tree-level mixing. 

We turn now to to the neutralino sector. When the lightest neutralino is bino- or 
Higgsino-like in this model, it shares the common features of the analogous neutralino 
in the MSSM. However, the masses differ at the SUSY scale for the same GUT 
scale parameters mainly as a consequence of gauge kinetic mixing in the gaugino 
sector which in our example amounts to shift of about 10 per-cent. This effect is 
especially important in the case of a bino LSP because the off-diagonal gaugino mass 
parameter M BB > can easily reach values of 0(— 50 GeV) (for Mi/ 2 ^ 200 GeV) or 
even 0(— 150 GeV) (for Mi/ 2 ^ 1000 GeV). For this purpose, we depict in the left 
column of Fig. 5 the mass of the lightest neutralino in the (m , iWi/^-plane as well 
the BLino fraction. To point out again the importance of kinetic mixing, we show the 
same plots in the right column of 5 without the kinetic mixing. A shift of masses is 
clearly visible and also the admixture differs by several orders. The other parameters 
are the same as for Fig. 4. 

For completeness we note that the differences between tree-level masses and loop- 
corrected masses are of similar size for the MSSM-particles. The R-sneutrinos receive 
somewhat larger loop corrections of about 3-4 per-cent compared to a few per-mille 
for the L-sneutrinos. However, interesting effects in the Higgs and neutralino sector 
can happen which we discuss in the next subsections. 

3.3 The Higgs sector 

In this section we concentrate on the Higgs sector and discuss new phenomenological 
aspects arising in the B — L model. In Fig. 6 we show the lightest scalar Higgs mass 
in the (m , Mi/ 2 )-plane. In addition, we give the bilepton fraction. The other input 
values are the same as for BLI. The nature of the second lightest Higgs is roughly 
complementary to the lightest one with respect to the ratio of the bilepton nature 
versus the Higgs doublet nature. The reason is that the mixing between the light 
states with the two heavy states, which have masses above a TeV, is quite small. 
Note that complete region shown is compatible with recent LHC data even though 
the mass of the second lightest Higgs boson is in most parts above 140 GeV, as it is 
mainly a bilepton with a small production cross section at the LHC. 

Close to the border of the allowed regions in the (m , Mi/ 2 )-plane, the lightest 
Higgs particles become bilepton-like. Not only can this be observed for a variation 
of mo and Mi/ 2 but also by adjusting tan/3', as shown in Fig. 7 where we have fixed 
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Figure 5: The (itiq, Mi/ 2 ) plane for a bino-LSP with kinetic mixing (left) and without 
(right). First row: mass of lightest neutralino. Second row: BLino fraction. The other 
parameters correspond to point BLI. 

m = 1000 GeV and My 2 = 500 GeV. To understand this behavior we neglect gauge 
kinetic mixing for simplicity as then the bilepton sector decouples from the MSSM 
Higgs bosons. In this limit we obtain at tree level a similar formula for masses as for 
the MSSM Higgs bosons: 

m h = 2 ( m z> + m Ao T ^{m 2 z , + m 2 A0 J 2 - ±m\,m\* cos 2 (2/3')) (3.1) 
Equations. (2.23) and (2.34) imply that for fixed Mz', rn 2 A0 shows a sizable depen- 

rj 

dence on tan /3'. We checked that very light bilepton-like Higgs scalars are not ruled 
out by experimental data using HiggsBounds 3.6.1beta [45, 46]. However, the 
mixing between the bilepton and the MSSM-like Higgs is rather small and thus the 
branching ratio h 2 — » h\h\ is at most a few per-cent. Therefore, the main decay 
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Figure 6: Mass of the lightest Higgs (first row on the left) and the second lightest Higgs 
(first row on the right), the logarithm of the bilepton fraction of the lightest Higgs (second 
row on the left) as well as the bounds from Higgs searches (second row on the right) 
in the (mo, Mi / 2 )-plane. For the Higgs search the saturation of the tightest bound as 
calculated by HiggsBounds is shown. While the blue is still consistent with all data, the 
red area is excluded. The most stringent channels are e + e~ — > Zh\ t 2,hx t 2 — > bb and 
pp — > h± t 2 — > W + W~ . The other parameters for point BLI. 

channels of the doublet Higgs are still SM final states and the well-known bounds 
do hold. However, as can be seen also in Fig. 7, the mass of the MSSM-like Higgs 
bosons gets pushed to larger values for very light bilepton scalar s. Such a behavior 
has already been observed in the literature when considering models with extended 
gauge symmetries [47, 48, 49, 50, 51, 52]. 

In Fig. 8 we take the point BLII and vary m and Mi/ 2 . We see that there is 
a sizable region where the lightest Higgs, being essentially a bilepton, has a mass 
of less than half of the second lightest, which is mainly like the MSSM h°. Even 
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Figure 7: a) masses of two lightest scalars. b) doublet (dashed blue) and bilepton (black) 
fraction of lightest higgs as function of tan/3'. The other input parameters are as for point 
BLII, but with m = 2M 1/2 = 1 TeV 



though the bilepton has only a small admixture of the doublet Higgs bosons, it is 
large enough to determine its main decay properties, which are mainly SM-like with 
respect to its decay into SM fermions. 

Loop corrections Concerning the loop-corrections in the Higgs sector, the picture 
is often comparable with the MSSM: the lightest two Higgs bosons receive very large 
corrections. These are even larger for the chosen BLSSM point than in the MSSM. 
For the heavy doublet scalar as well as the charged and MSSM-like pseudoscalar 
Higgs, the differences between tree-level and one-loop masses are rather small and of 
the same size as in the MSSM. The neutrino Yukawa couplings don't play any role 
in this context because the correct explanation of neutrino data requires them to be 
very small. However, there are sizable loop corrections due to the large Y x couplings 
similar to the top-stop contributions to the lighter MSSM-Higgs. For example, in 
the case of the BLIV scenario, m^ 2 gets shifted from about 252 GeV at tree level to 
about 210 GeV at the one-loop level. This is a consequence of the mass hierarchy 
between fermions and bosons in the extended gauge/Higgs sector. 

For completeness we note that in the Higgs sector one finds that the mass of the 
bilepton-like Higgs field vanishes at tree level in the limit tan /3' — > 1. Note, however, 
that tan /3' = 1 is a saddle point and not a minimum of the tree-level potential. This 
in turn implies that the loop corrections will be large compared to the tree level 
similar to as it is in the MSSM when considering there the limit tan /3 — > 1. We 
explicitly demonstrate this in Fig. 9 where we compare the tree-level and one-loop 
masses of the two lightest Higgs fields, fixing the input parameters as in point BLII 
but varying tan/3'. This behavior is also reflected in the values of fi calculated from 
the one-loop and tree-level tadpole equations as shown in the right plot of Fig. 9. 
It can be understood from eq. (2.20): compared to cos(2/3'), all parameters and the 
one-loop corrections show only a very mild dependence on /3' as tan /3' has to be close 
to one. Therefore, the one-loop correction to the tadpole equation can be included 
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Figure 8: Mass of the two lightest Higgs fields (first row) as well as the logarithm of the 
bilepton fraction (left plot in second row) in the (mo, Mi/2)-plane. The right plot in the 
second row shows the saturation of the tightest bound as calculated by HiggsBounds: the 
blue are is allowed, the red one excluded by Higgs searches. The most sensitive channels 
are e + e~ — > ZI12, 112 — > bb, pp — > A — > tt and pp — > /12 — > W + W~ . The other parameters 
are based on BLII. 

in the first term and effectively be absorbed in a redefmtion of /3', denoted by (3', 
using the equation ggBL% 2 cos (2/3') + 4 ^.~^ = ggsLX 2 cos(2^') with a shifted (3'. 

In Fig. 10 we show the masses of the two lightest Higgs fields for BLI with a 
variation of Mz>- As expected the tree-level and one-loop mass of the lightest Higgs 
which is consists mainly of the SU(2) doublet is nearly independent of Mz>- In 
contrast to that, the tree-level mass of the bilepton-like Higgs depends strongly on 
Mz>- Furthermore, the one-loop corrections can be nearly of the order known from 
the MSSM for the doublet Higgs depending on the mass of the Z'. Also the sign of 
the correction can change depending on the mass ordering of Z' and the BLino-like 
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Figure 9: Left: mass of the lightest Higgs fields at one- loop level (black) and tree level 
(dashed blue). Right: value of /u calculated from the one- loop corrected tadpole equations 
(black) and the tadpole equations at tree level (dashed blue). The other input parameters 
correspond to those of Fig. 7 
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Figure 10: Left: masses of the lightest Higgs (dashed) and second lightest Higgs (solid) 
at one-loop level (black) and tree level (blue) for a variation of Mz>- Right: relative 
difference between tree-level and one-loop mass for the lightest Higgs (dashed blue) and 
second lightest Higgs (solid black). The other input parameters correspond to point BLI. 



neutralino. 



3.4 The neutralino sector 

Similarly to how it is in the CMSSM, the lightest neutralino is often bino-like and 
the main difference is, in this case, that the relation between the parameters at 
different scales gets changed due to the gauge kinetic mixing. Note that this holds 
even though the soft-breaking gaugino mass term Mb> is always smaller than Mi, 
because, at one-loop level and without kinetic mixing, the relation 

Mi/2 _M 1 _M K . , 

2 — 2 — 2 \~-^) 

9gut 9y 9bl 

would hold and gsh is always smaller than gy if unification at the GUT scale is 
assumed, as can be seen in eq. (2.7). However, usually there is a large mixing 
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Figure 11: The content of the LSP depending on the tan/3' (right) and Mz> (left). 
The other parameters are those of BLIII. The color code is as follows: gaugino fraction 
(dashed black), Higgsino fraction (blue), logarithm of the BLino fraction (dotted black) 
and bileptino fraction (dot-dashed blue). 

between the BLino with the bileptinos, leading to heavy states. However, we will 
demonstrate that nevertheless regions exist where the lightest neutralino is BLino- 
or even bileptmo-like. Therefore, a neutralino LSP can have four different natures 
in the BLSSM (bino, Higgsino, BLino, bileptino) in contrast to only two possibilities 
in the CMSSM. This can provide interesting features in the context of dark matter 
[34]. 

It is well known that in the MSSM, it is very hard to reach Higgsino fractions of 
the LSP larger than 50 per-cent: even in the focus point region this fraction hardly 
ever exceeds 30 per-cent. Only in a tiny region where |//| gets close to 0, and which is 
excluded by LEP data, does it get larger than 50 per-cent. In contrast in our model 
new contributions show up in the formula for /i in eq. (2.20), in particular the term 
ggBLX 2 cos(2/3'). Using universal boundary conditions it is in general negative and is 
particularly sizable for large tan /3' and M^/. Therefore it is possible to increase the 
Higgsino fraction of the LSP by increasing these two parameters with fixed CMSSM 
parameters as shown in Fig. 11. 

As mentioned above, the soft-breaking parameter Mb' is always smaller than 
Mi, but the large mixing between the BLino and the bileptino usually implies that 
the bino is still the LSP. However, there are regions where this mixing is small and 
the BLino becomes the LSP. In particular this happens if // 3> Qbl^ — Mz> which 
happens either for large \Y X \ or large m , as this increases the difference m? — m 2 . As 
an example we show in Fig. 12 the dependence of //, the masses of all neutralinos and 
the content of the lightest on mo. As claimed, // grows with increasing mo leading 
to a larger mass splitting between the bileptino-like neutralinos and the others. For 
very large values of //, the bilepton fields are nearly decoupled and the nature of the 
LSP becomes BLino-like. In this case one has to check if one can obtain the correct 
value for the relic density. A principal possibility are resonances as there are two 
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Figure 12: a) fjf as function of tuq. b) masses of all neutralinos. c) content of the 
lightest neutralino: gaugino fraction (dashed black), Higgsino fraction (blue), logarithm 
of the BLino fraction (dotted black) and bileptino fraction (dot-dashed blue). The input 
parameters were those of BLIV but with M 1 / 2 = 1 TeV. 

light Higgs bosons and the LSP mass could easily be half of one of the Higgs masses. 
However, it still has to be checked if the corresponding couplings are sufficiently 
large, which however is beyond the scope of this paper. 

We want to close the discussion of the BLino LSP with a remark about the 
importance of the loop corrections. It is well known that in the MSSM one gets a 
few per-cent corrections to the masses of sleptons, neutralinos and charginos [26]. 
In the model considered, the corrections are usually of a similar size. However, this 
doesn't apply for a light BLino because here loops contribute with rather heavy 
particles, in particular A® and the bileptino-like neutralinos. This is demonstrated 
by inspecting scenario BLIV. Varying tan /3' we find that these corrections get larger 
the closer tan/3' gets to one as can be seen in Fig. 13. At tree level one might 
then conclude that the LSP could be massless even for unified gaugino masses. This 
behavior can be roughly understood when neglecting gauge kinetic mixing as then 
the bileptinos and the BLino decouple from the MSSM neutralinos. In this limit it 
is thus sufficient to consider only the lower left 3x3 block of the neutralino mass 
matrix given in eq. (2.40). Taking tan/3' 1 or equivalently v v ~ Vfj we find for the 
LSP mass 

1 -{M B + // - ^(M B - fx') 2 + 8(g BL v v ) 2 ) (3.3) 
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tan /?' tan /?' 

Figure 13: Left: mass of lightest neutralino at tree level (dashed blue) and one-loop level 

T 

(black) for a variation of tan/3'. Right: relative size of the correction |Am| = |1 — I- 
The input parameters are those of BLIV. 



This expression can obviously be negative or positive for the same sign of // depend- 
ing on the value of v v . Although this gets changed at the one-loop level, it is still 
easy to obtain a dark matter candidate within the mass range preferred by direct 
detection experiments like DAMA [53, 54], CRESST [55] or Cogent [56]. 
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Figure 14: LSP with large bileptino fraction (benchmark scenario BLV): a) mass of 
neutralinos, b) neutralino content. The color code on the right-hand side is as follows: 
gaugino fraction (dashed black), Higgsino fraction (blue), logarithm of the BLino fraction 
(dotted black), bileptino fraction (dot-dashed blue). 

Finally, we note that also a bileptino-like LSP can be obtained in this model. 
The necessary condition, \fi'\ being smaller than and all gaugino mass parameters, 
can be obtained if the difference between and m? becomes small. This can be 
accommodated by adjusting the entries of Y x . As an example, we show in Fig. 14 
the masses of all neutralinos as well as the composition of the lightest neutralino 
as function of Y Xj u while keeping all other values as in scenario BLV. Already a 
10 per-cent decrease leads to a nearly a pure bileptino LSP and its mass depends 
strongly on Y x u. For larger values a level crossing takes place and the LSP becomes 
bino-like. In principal this coupling could be larger at the electroweak scale but then 
one would encounter a Landau pole below the GUT scale. 
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4. Conclusions and discussion 

We have discussed in this paper the mass spectrum of the minimal B — L extension 
of the MSSM taking universal boundary conditions at the GUT scale. We have cal- 
culated the spectrum using two-loop RGEs and the complete one-loop contributions 
to all masses, which are particularly important in the Higgs and neutralino sectors. 
Consistency with current bounds on the additional Z' implies that the scalar part- 
ners of the fermions are quite heavy in this scenario, except the sneutrinos can be 
light under certain conditions. However, this is a consequence of the taking the mass 
parameters of the sfermions equal to the ones in the extended Higgs sector. Relax- 
ing this assumption allows for lighter non-sneutrino sfermions in addition to light 
sneutrinos. 

It turns out that gauge kinetic mixing between the two Abelian gauge groups is 
quite important for Higgs bosons and neutralinos and it cannot be neglected. On 
one hand it leads to sizable shifts in the masses of up to 10 per-cent. On the other it 
induces tree-level mixing between the MSSM states and the states of the extended 
B — L sector leading to important shifts in the nature of the corresponding particles. 
This holds in particular for the light Higgs bosons and the lightest neutralino. For 
example in the latter case we find regions in parameter space with light neutralinos 
as preferred by DAMA or COGENT. Moreover, the nature of the lightest neutralino 
can be quite different from the usual CMSSM, e.g. we have identified regions where 
it is either dominantly Higgsino-, BLino- or bileptino-like. 

In the extended Higgs sector we find that one-loop corrections are not only 
important for the MSSM-like h° but also for the light bilepton field. This particle 
can be so light that the MSSM ft°-like state can decay into two of them without 
conflicting with any of the known experimental results. However, in general we find 
that the corresponding branching ratio is at most a few per-cent. 
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A. Mass matrices 

Here we collect the tree-level formulas for the remaining sfermion mass matrices. 
• Mass matrix for Sleptons, Basis: {e L ,e R ) 




(A.l) 
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mLL = m 2 L + ^Y e + i [[g\ - g\ + g 2 + gg BL )(v 2 d - v 2 u ) + 2(gg BL + g 2 BL )(v 2 - 

(A.2) 

man = m 2 E + ^Y e Y} + l - {{2g 2 + 2 ^ 2 + ggsh) (v 2 u - v 2 d ) - 2 (2^bl + «^ L ) (t; 2 



(A.3) 



Mass matrix for Down-Squarks, Basis: (^d L ,d,Rj 



/ "' / 1 ^(^J-^yJ) m M 1 



m LL = + ^y d ty d + ^ ((s 2 + 3g 2 + ~g 2 + ~gg BL ){v 2 u - v 2 d ) + 2(g 2 BL + ^ B l)(^ - <)) 3 

(A.5) 

m** = m 2 D + + 1 ( - (2<? 2 + 2g 2 - gg BL )(vj - v 2 u ) + 2{g 2 BL - 2gg BL )(v 2 - v 2 )l 



(A.6) 



Mass matrix for Up-Squarks, Basis: (ul,u b ) 



m|=( m M Uv,T,-v^)^ (A7) 



m LL = m 2 Q + + — ([g\ - 3g 2 2 + ~g 2 + ^ L )(t; 2 - ^) + 2(^ Bi + g 2 BL )(v 2 - v 2 )) 

(A.8) 

m RR = m 2 v + V -^Y u Yl + — (2(g 2 BL + A~gg BL )(v 2 - v 2 ) + {Agl + Ag 2 + gg BL ){v 2 d - v 2 u )) 1 

(A.9) 

B. RGEs 

The calculation of the renormalization group equations performed by SARAH is based 
on the generic expression of [40]. In addition, the results of [22] are used to include 
the effect of kinetic mixing. 

The (3 functions for the parameters of a general superpotential written as 

W(<f>) = ^r'rw, + \.Y l!l O,0 J 0, (B.l) 
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can be easily obtained from the shown results for the anomalous dimensions by using 
the relations [57, 58] 



(3f = Y p{ - ij lp V , (B.2) 
= ^7/ . (B.3) 

For the results of the other parameters as well as for the two-loop results which we 
skip here because of their length we suggest to use the function CalcRGEs [] of SARAH 
with the model files shown in appendix C. 

B.l Anomalous dimensions 



- ^ (2(<?yy + 9yb) + 2VlO(g YY g BY + QybQbb) + 5 (lfyf + 32gj + g 2 BB + g 2 BY ^ ) 1 

(B.4) 

7 f ) = yX + y;y v t 

- ^ (2(^yy + g YB ) + 2\f\Q(g Y YgBY + gY BgBB) + 5 (2g\ + g 2 BB + g BY ^j ) l) 

(B.5) 

7 g = 3Tr (Y d Yj) + Tr (y e Y^ - A (bgj + g 2 YY + g 2 YB ) (B.6) 
7 g = 3Tr (YX) + Tr (Vj) ~ ^ (5s 2 2 + <? yy + 5yB ) (B.7) 

+ A (W 10(# yy # By + gY bQbb) - 5 (32# 2 + + g 2 BY ^j _ 8(# 2 y + g 2 YB ) ) 1 

(B.8) 

(1) = 2r *yT 

--(32(s yy + £ yB ) + 5(32^ 2 + 

fl'ss + fi'ly) + 8v / 10(fi'Yyfi'By + gY b9bb)^ 1 

(B.9) 

7? } = 2F;F e T - A ^(gyy^y + M b) + 5 (<?| B + <?| y ) + 8(g 2 YY + <? yB )) 1 



(B.10) 

7^ = 2YX + 2{YX + Y X *Y X ) - ^ («£ B + g 2 BY y (B.ll) 
7? } = -3(sIb + <?|y) + 2Tr(y^t) (B.12) 
7i 1} = -3(^ + ^y) (B.13) 
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B.2 Gauge Couplings 

We give here and in the subsequent section the beta functions for the RGEs of the 
gauge couplings and gaugino mass parameters in a basis independent way as the 
off-diagonal values for the U(l) gauge couplings and gaugino mass paramters are 
generated due to RGE effects. 
3 



= - (llg YY + 4Vl0g YY g BY + g YY (llg YB + 15g 2 BY + 2V^9yb9bb) 
+ g YB (l5g BB + 2Vldg YB ^g BY ^ (B.14) 

Psibb = 5 [} 1 9yb9bb + AV\Qg YB g BB + 15^| B + Wg YY g YB g BY + 2y/l0g YY g BB g BY 
+ 2VWg YB g 2 BY + l5g BB g 2 BY ) (B. 15) 



yb = 5 {pyy (i59bb + 2VWg YB ^)g BY + g YY (llg YB + 2Vl6g BB ^j 

+ g YB (llg 2 YB + 15g% B + 4Vl0g YB g BB )) (B.16) 



PgsY = 5 ( U 9yy9by + 9yy{i19yb9bb + 2Vw(2g 2 BY + g 2 BB ^ 

+ g BY (l5(g 2 BB + g 2 BY ^+2VWg YB g BB ^ (B.17) 



P$?=92 (B.18) 
P£ = -3^ (B.19) 



B.3 Gaugino Mass Parameters 

(3^1 = - (llgYyMi + g BY (l5g BB M BB/ + lhg BY M 1 + 2Vl0g YB M BB ?) 

+ g YY (llg YB M BB , + 2V1(WMbb' + 4Vl(WMi)) (B.20) 

ft™ = 2g\M 2 (B.21) 

0® = -6g 2 3 M 3 (B.22) 

^ = -(llg 2 YB M B + 15^bb(^bbM b + g BY M BB ^ + 

2VWg YB (2g BB M B + g BY M BB ^j + # yy (ll#y B + 2^^) M BB ^ (B.23) 
= — ^11 g YY M BB / + llg YB M BB / + 2y/lOg YB (2g BB M BB > + g BY (m x + M b )) + 
15 {g 2 BB M BB , + g BB g BY (m x + M B ) + g 2 BY M BB ^ 
+ 9yy (ll^B (Mi + M B ) + 2v / ld(2^yM B B' + 9bb (m x + M 
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C. Model files for SARAH 



Below we list the model files used for SARAH to study the model presented in this 
paper. Using this one can generate the Fortran code for the SPheno extension to 
reproduce the results presented. These files will also become part of the public 
SARAH package in near future. 

C.l B-L-SSM.m 



(* 

(* Particle Content*) 

(* 

(* Gauge Superfields *) 

Gauge [[ 1 ]] = { B, U[l], hypercharge , gl, False}; 

Gaugc[[2]] = {WB, SU[2], left, g2 , True } ; 

Gaugc[[3]] = {G, SU[3] , color, g3, False}; 

Gauge [ [4] ] = { Bp , U[l] , BminusL , glp , False}; 

(* Chiral Superfields *) 



Fields 


[1]] 


= {{uL, dL}, 


3, q, 1/6, 2, 


3, 1/6}; 


Fields 


[2]] 


= {{vL, cL}, 


3,1, -1/2, 2, 


l, -1/2}; 


Fields 


[3]] 


= {{HdO, Hdm}, 


1, Hd, -1/2, 2, 


1, 0}; 


Fields 


[4]] 


= {{Hup, HuO}, 


1, Hu, 1/2, 2, 


1, 0}; 


Fields 


[5]] 


= { conj [dR] , 3 , 


d, 1/3, 1, -3 


-1/6}; 


Fields 


[6]] 


= {conj [uR] , 3 , 


u, -2/3, 1, -3 


-1/6}; 


Fields 


[7]] 


= {conj [cR] , 3 , 


e , 1,1, 1 


1/2}; 


Fields 


[8]] 


= {conj [vR] , 3 , 


vR, 0, 1, 1 


1/2}; 


Fields 


[9]] 


= {CIO, 1, CI, 


0, 1, 1, -1}; 




Fields 


[10] 


] = {C20, 1, C2, 


0, 1, 1, 1}; 





■*) 
■*) 



(* *) 

(* Superpotential *) 

(* *) 

SupcrPotential = { {{1, Yu} , { u , q , Hu} } , { { - 1 ,Yd} , { d , q , Hd} } , 

{{-l,Yc},{e,l ,Hd}}, {{l,\[Mu]},{Hu,Hd}}, 

{{l,Yv},{l ,Hu,vR}}, {{-l,MuP},{Cl,C2}}, 
{{l,Yn},{vR,Cl,vR}} }; 

(* *) 

(* Integrate Out or Delete Particles *) 

(* *) 

IntegrateOut={}; 
DeleteParticles={}; 

(* *) 

(* ROTATIONS *) 

(* *) 

(* Different eigenstates: gauge eigenstates and eigenstates after EWSB 

*) 



32 



NameOfStates={GaugcES , EWSB} ; 

(* Gauge fixing terms for Gauge eigenstates 

DEFINITION [ GaugcES ] [ GaugcFixing] = 

{ {Dcr[VWB], -1/(2 RXi [W] ) } , 
{Dcr[VG] , -1/(2 RXi [G] ) }}; 



*) 



(*" 



Rotations in gauge sector 



*) 



DEFINITION [EWSB] [ GaugeSector ] = 
{ {{VB,VWB[3] ,VBp} ,{VP,VZ,VZp} ,ZZ} , 
{{VWB[1] ,VWB[2] } ,{VWm, conj [\Wm] } ,ZW} , 
{{fWB[l] ,fWB[2] ,fWB[3]} ,{fWm,fWp,fWO} ,ZfW}}; 



VEVs 



DEFINITION [EWSB] [VEVs] = 
{{SHdO, {vd, 1/Sqrt[2]} 
{SHuO, {vu, 1/Sqrt[2]} 
{SvL, {0, 0}, {sigmaL, 
{SvR, {0, 0}, {sigmaR, 
{SC10, {xl, 1/Sqrt[2]} 
{SC20, {x2, 1/Sqrt[2]} 



(* Matter Sector 



{ sigmad , \[ImaginaryI]/Sqrt[2]} ,{phid,l/Sqrt[2]}} , 
{sigmau , \[ImaginaryI]/Sqrt[2]} ,{phiu,l/Sqrt[2]}} , 
\[ImaginaryI]/ Sqrt [2] } , { phiL , 1/ Sqrt [2] } } , 
\[ImaginaryI]/Sqrt [2]} , { phiR , 1 / Sqrt [2]}} , 
{sigmal, \[ImaginaryI]/Sqrt [2]} ,{phil , 1/Sqrt[2]}}, 
{sigma2, \[ImaginaryI]/Sqrt [2]} ,{phi2 , 1/ Sqrt [ 2 ] } } } ; 



*) 



DEFINITION [EWSB] [ MattcrScctor] = 
{ {{SdL, SdR}, {Sd, ZD}}, 

{{SuL, SuR}, {Su, ZU}}, 

{{ScL, ScR}, {Sc, ZE}}, 

{{sigmaL , sigmaR} , {Svlm, ZVI}} , 

{{phiL ,phiR} , {SvRe, ZVR} } , 

{{phid, phiu,phil, phi2}, {hh, ZH}}, 

{{sigmad, sigmau , sigmal , sigma2}, {Ah, ZA}}, 

{{SHdm, conj [SHup] } ,{Hpm,ZP}} , 

{{fB, fWO, FHdO , FHuO.fBp ,FC10,FC20} , {L0, ZN}}, 
{{{fWm, FHdm}, {fWp, FHup}}, {{Lm,UM} , {Lp,UP}}}, 
{{FvL.conj [FvR] } , {Fvm,UV} } , 

{{{FcL} ,{conj [FcR] } } ,{{FEL,ZEL} ,{FER,ZER}}} , 
{ { { FdL } ,{ conj [FdR] }} ,{{FDL,ZDL} ,{FDR,ZDR} } } , 
{{{FuL} ,{conj [FuR]}} ,{{FUL,ZUL} ,{FUR,ZUR}}} 

}; 



Gauge Fixing after EWSB 







DEFINITION [EWSB] [ GaugeFixing] = 
{ {Dcr [VP] , 

{Der[VWhi]+\[ImaginaryI] Mass [Win] RXi [W] Hpm[{l}] 
{Dcr[VZ] - Mass[VZ] RXi[Z] Ah[{l}], 
{Der[VZp] - Mass [VZp] RXi[Zp] Ah[{2}], 
{Dcr [VG] , 



1/(2 RXi[P])}, 
l/(RXi[W])} , 
1/(2 RXi[Z])}, 
1/(2 RXi[Zp])}, 
1/(2 RXi [G] ) } } ; 



Phases 



*) 
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DEFINITION [EWSB] [ Phascs] = 
{ {fG, PhaseGlu} }; 



-*) 
*) 
-*) 

(* Dirac Spinors for gauge eigenstates *) 

DEFINITION [ GaugcES ] [ DiracSpinors]={ 
Bino ->{fB , conj [ f B ] } , 
Wino -> {fWB, conj [fWB] } , 
Glu -> {fG, conj [fG]} , 
HO -> {FHdO , conj [FHuO]} , 
HC -> {FHdm, conj[FHup]}, 
Fdl -> {FdL , 0} , 
Fd2 -> {0, FdR} , 
Ful -> {FuL, 0} , 
Fu2 -> {0 , FuR} , 
Fel -> {FeL, 0} , 
Fe2 -> {0 , FeR} , 
Fvl -> {FvL , 0} , 
Fv2 -> {0 , FvR} , 
FC -> {FC10, conj[FC20]}, 
FB -> {fBp , conj [fBp] } 

}; 

(* Dirac Spinors for eigenstates after EWSB *) 

DEFINITION [EWSB] [ DiracSpinors] = { 
Fd ->{ FDL, conj [FDR]}, 
Fc ->{ FEL, conj [FER]}, 
Fu ->{ FUL, conj [FUR]} , 
Fv — >{ Fvm, conj [Fvm]} , 
Chi ->{ L0 , conj [L0] } , 
Cha — >{ Lm, conj [Lp]} , 
Glu ->{ fG, conj [fG]} 

}; 



(* 

(* Dirac Spinors 

(* 



C.2 SPheno.m 



(* *) 

( * MINPAR * ) 

(* *) 

MTNPAR={{l,mO} , 
{2,ml2}, 
{3 .TanBcta} , 
{4 ,SignumMu} , 
{ 5 , Azcro } , 
{ 6 , SignumMuP } , 
{7 ,TanBetaP} , 
{8,MZp}}; 

RealParameters = {TanBeta, TanBetaP } ; 



Tadpoles and renormalization scale *) 
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(* *) 

ParametcrsToSolveTadpoles = {B[\[Mu]] ,B[MuP] ,\[Mu] ,MuP}; 

RenormalizationScaleFirstGuess = m0"2 + 4 ml2~2; 
RenormalizationScale = MSu [ 1 ] *MSu [ 6 ] ; 



(* 

(* Boundary conditions 

% (* 



(* Definition of GUT scale *) 

ConditionGUTscale = ( gl * glp-glglp* glpgl ) / Sqrt [ glp '2+glpgl " 2] = g2 

(* Boundary conditions at GUT scale *) 

Boundary HighScale={ 

{gl ,(gl*glp-glglp*glpgl)/Sqrt [glp "2- glpgl "2]} , 

{gl,Sqrt[(gr2+g2-2)/2]}, 

{g2,gl}, 

{glp, gl}, 

{glglp,0} , 

{glpgl ,0} , 

{T[Ye] , Azero*Yc} , 

{T [ Yd] , Azcro*Yd} , 

{T[Yu] , Azero*Yu} , 

{T[Yv] , Azcro*Yv} , 

{T[Yn] , Azcro*Yn} , 

{mq2 , DIAGONAL mO " 2 } , 

{ ml2 , DIAGONAL mO " 2 } , 

{md2 , DIAGONAL mO " 2 } , 

{mu2 , DIAGONAL mO " 2 } , 

{ me2 , DIAGONAL mO ' 2} , 

{mvR2 , DIAGONAL mO " 2 } , 

{mHd2 , mO " 2 } , 

{mHu2 , mO " 2 } , 

{mC12, m0"2}, 

{mC22, m0"2}, 

{MassB, ml2}, 

{MassWB,ml2} , 

{MassG,ml2} , 

{MassBp ,ml2} , 

{MassBBp.O} , 

{MassBpB,0} 

}; 



(* Boundary conditions at SUSY scale *) 

BoundarySUSYScalc = { 

{glT,(gl*glp-glglp*glpgl)/Sqrt [glp'2+glpgl "2]} , 

{glpT, Sqrt[glp"2+glpgl-2]}, 

{glglpT,(glglp*glp+glpgl*gl)/Sqrt [glp-2+glpgl -2]} , 

{gl, glT}, 

{glp, glpT}, 

{glglp, glglpT}, 

{glpgl ,0} , 

{vcvP, MZp/glp}, 

{ bctaP , ArcTan [ TanBetaP ] } , 

{x2 , vevP*Cos [betaP] } , 

{xl , vevP*Sin [betaP] } , 

{Yv, LHInput[Yv]} , 

{Yn, LHInput[Yn]} 
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}; 

(* Boundary conditions at EWSB scale *) 

BoundaryEWSBScale = { 

{glT,(gl*glp-glglp*glpgl)/Sqrt [glp'2+glpgl "2]} , 

{glpT, Sqrt[glp"2+glpgl-2]}, 

{glglpT ,(glglp*glp+glpgl*gl)/Sqrt [glp-2+glpgl -2]} , 

{gl, glT}, 

{glp, glpT}, 

{glglp, glglpT}, 

{glpgl ,0} , 

{vevP, MZp/glp}, 

{ bctaP , ArcTan [ TanBctaP ] } , 

{x2 , vcvP*Cos [betaP] } , 

{xl , vcvP*Sin [bctaP] } 

}; 

(* Initialization values *) 

InitializationValues = { 

{glp, 0.5}, 

{glglp, -0.06}, 

{glpgl, -0.06} 

} 

(* Boundary conditions for SUSY scale input *) 

Boundary LowScaleInput={ 
{vd , Sqrt [4 mz2/ ( gl ~2+g2 " 2)] * Cos [ ArcTan [ TanBcta ] ] } , 
{vu , Sqrt [4 mz2/(gl"2+g2 " 2)] * Sin [ ArcTan [ TanBeta ] ] } 

}; 

(* *) 

(* Two and Three body decays *) 
(* *) 

ListDecayParticles = Automatic; 
ListDecayParticlcs3B =Automatic ; 
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Table 2: Points with distinct features: BLI is similar to the MSSM with a bino LSP, BLII 
provides a very light, bilepton-like Higgs boson, BLIII has a very large Higgsino fraction, 
BLIV has a BLino LSP and sizable mixing between the doublets and bileptons in the Higgs 
sector, and BLV has a bileptino LSP. h$M denotes the Higgs fields which is most similar 
to the SM Higgs particle. 
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Table 3: Comparison of mass spectrum in the MSSM and the BLSSM for input parameters 
BLI of Tab. 2. We give the masses for one- and two-loop RGE-evaluations. In addition, 
we include a comparison of the case of properly taking into account gauge kinetic mixing 
(KM) versus neglecting it (NKM). 
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